Impedance of the single electron transistor at radio-frequencies 
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We experimentally characterise the impedance of a single electron transistor (SET) 
at an excitation frequency comparable to the electron tunnel rate. Differently from 
usual rf-SET operations, the excitation signal is applied to the gate of the device. At 
zero source-drain bias the single electron transistor displays both resistive (Sisyphus 
resistance) and reactive (tunnelling capacitance) components to its impedance. We 
study the bias dependence of the complex impedance, investigating its response as 
the electron tunnel rate becomes large with respect to the driving frequency. The 
experimental data are compared to values calculated from a master equation model. 

PACS numbers: Valid PACS appear here 

Keywords: single electron transistor, rf-SET, Sisyphus resistance 



a ) cc538@cam.ac.uk 

b ) ajfl006@cam.ac.uk 



1 



The single electron transistor (SET) is a fundamental nano-scale electronic device, it con- 
sists of an island coupled to source and drain leads by two low-capacitance tunnel junctions.- 
Its conductance is modulated by the polarization charge induced on the island by its electro- 
static environment. As a result, it is possible to use the SET as an ultra-sensitive and high- 
bandwidth charge transducer, capable of detecting single electrons on a sub-microsecond 
timescale.- - - An important characteristic that determines the performance of a SET is the 
rate with which it responds to an external perturbation. This is limited by the tunnelling 
rate of the electrons at the source and drain junctions and determines the ultimate charge 
detection bandwidth. When the time scale of the perturbation approaches the tunnelling 
time of the electrons, the SET response becomes susceptible to the stochastic nature of the 
tunnelling events. In this Letter we quantify the complex impedance that arises due to the 
competition between tunnelling rates and radio-frequency excitation. 

The dissipative response (Sisyphus resistance) of a single electron box (SEB}2 to a radio- 
frequency signal has previously been measured. A single electron device can also exhibit a 
capacitance due to electron tunnelling events leading, rather than being exactly in phase 
with, the excitation. This capacitance has been investigated in gallium arsenide quan- 
tum dots^ and, similar to the Sisyphus resistance, allows measurements of single electron 
charging with just a single tunnel junction. This tunnelling capacitance is distinct from 
the quantum capacitance that arises from bandstructure curvature in single Cooper pair 
devices^ 1 ^ and double quantum dots.— 

Aluminium SETs were fabricated by a standard double angle evaporation technique using 
bilayer resist and controlled oxidation. 13 They were measured in a dilution fridge at an 
electron temperature lower than 200 niK. A magnetic field (B=600 mT) was applied to 
suppress superconductivity in the aluminium. A radio frequency resonant circuit, consisting 
of a chip inductor, was connected to the gate of the SET (Fig. 1(a) & (b)). This is different 
from the usual configuration of the rf-SET 2 , where the resonant circuit is connected to the 
source-drain of the SET, and where the largest contribution to the reflection coefficient is 
from modulation of the differential conductance. 

The circuit was driven at resonance (/o = 347 MHz). The amplitude of the rf-signal 
sent to the device was dn g = C g dV g ~ 0.09 e (2xQx<5V^xC 5 =2x60x2/iKx0.1fF, where 5V g 
is the amplitude of the rf-drive and Q is the resonator's Q-factor, estimated from the 3dB 
part of the resonance curve shown in FigfjTc)). After amplification by a low temperature 
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FIG. 1. (Color online) (a) A schematic of the resonant circuit. The SET is connected 
from the gate to the chip inductor (L = 560 nH) which, together with a parasitic 
capacitance (C p = 0.4 pF), comprises the resonator. A bias-tee on the PCB allows 
us to apply a dc voltage to the gate. Chip capacitors (Co = 150 pF) are connected 
from the source and drain to ground to provide an rf ground, (b) Scanning electron 
micrograph of the aluminium SET. Sample resistance at low temperature is 2R = 
200/cS7, where R is the resistance of each tunneling junction, (c) Amplitude and 
phase characteristic of the reflected signal at V s d = QV as a function of the frequency. 

and room temperature amplifier, the reflected signal from the resonant circuit was mixed 
with the reference signal to provide phase sensitive detection. The output of the mixer was 
amplified, low pass filtered and digitised with an oscilloscope. By using a line stretcher we 
varied the phase of the incident signal to the device with respect to the reference signal. As 
a result, we obtain the in-phase and quadrature components of the reflected signal, allowing 
us to calculate the amplitude and phase response of the SET. 

We briefly compare the differential conductance of the SET, measured using a standard 
low-frequency lock-in technique, to the amplitude and phase response of the reflected signal 
(Fig. [2]). In the phase response we can clearly observe the Coulomb diamonds that appear 
in the differential conductance graph, these corresponding to the blockade regime in which 
transport through the island is forbidden. Additionally there is a peak in the phase shift 
around the zero bias degeneracy points {V s d — V, n g — (m + 1/2) e). In the amplitude 
response, there is a dissipative signal at the degeneracy points that reduces with increasing 
source drain bias. Strikingly, by comparison with the phase response, the Coulomb diamonds 
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FIG. 2. (Color online) (a) The differential conductance of our SET is measured using 
a low-frequency lock-in amplifier. The total capacitance of the island Cs = 0.6 fF, 
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corresponds to a charging energy of E c = = 0.13 meV, the gate capacitance is 
C g = 0.1 fF. (b) Phase shift of the reflected signal, (c) Amplitude of the reflected 
signal relative to the refection in the blockaded region. 

are not observed. 

We now analyse the phase response, starting with observation of the Coulomb diamonds. 
When the device is in Coulomb blockade (CB), the gate capacitance (C g ) and the junction 
capacitances {C\ + C 2 ) are in series (Figj3](a)), with little contribution to the impedance 
from the tunnel junction resistance. Outside the blockaded region, the capacitance becomes 
shunted by the junction resistance. This results in an increase in capacitance ( ^ 9 ^ 1 ~ 1 ^^ ) 
seen by the resonator, lowering the resonant frequency, and explaining the phase shift. 

Near the zero bias degeneracy points a description of the electron tunnelling dynamics is 
necessary. With the rf signal superimposed on the gate, the island cyclically transits between 
the nearly degenerate n and n + 1 charge states. If the tunnel rates at zero instantaneous rf 
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FIG. 3. (Color online) (a) The top left diagram shows the equivalent circuit for the 
two tunneling junctions in parallel. This includes the resistance and the geometrical 
capacitance of the junction. The top right image is a snapshot of the graph of the 
reflected signal phase as a function of the dc voltage applied to the gate and of the 
source to drain bias. On the graph three points are shown: A indicates the point 
of zero source to drain bias, B marks the peak in the phase shift at low source to 
drain biases and C indicates the region of high source to drain biases. The three 
diagrams at the bottom show the tunnelling events on the time scale of the rf-drive 
for each of these three points. I ne t is the net current to ground that results from 
the tunnelling of the electrons. We see that, as V s d is increased (B), the tunnelling 
events happen "early" in the rf-cycle and I ne t develops a negative phase shift with 
respect to the drive. If V s d is further increased (C), however, I net is generated by 
stochastic tunnelling events no longer correlated to the drive, (b) and (c) Graph of 
the measured and calculated variation of the reflected signal phase shift as a function 
of the dc gate voltage (b) and of the source to drain voltage (c). 



bias (we estimate at the zero bias degeneracies Y = ~2 GHz, where R is the resistance 
of one tunnelling junction), are greater than the drive frequency, the electrons on average 
tunnel early in the rf cycle. This phase shifts the current with respect to the rf voltage, 
leading to an effective capacitance C e ff in parallel with the SET. The total capacitance C 



can then be written as C = 1 c ^ 9 + ^+ Jy > , where < en > is the average charge on the 
island. This expression can be found from the expression of the polarisation charge induced 
on the island by the gate Q g = en — [C\(V s d — Vi) — C^Vi], where Vi is the island potential. 
The first term represents the dc limit of the capacitance, obtained by considering C g in 
series with the capacitances of the source and drain tunnelling junctions. The second term 
represents the contribution from C e ff and can be written as ~^ ^p P ^ 9 > where (1/C g )dn g 
is the amplitude of the rf drive and P n is the probability of having n electrons on the island. 
When k B T « E c and dn g « e, tunnelling to higher energy states can be neglected and 
P n is found by solving the Master equation which involves the states n and n + 1 1 



Pn — ^n+l,nPn+l ~ ~^n,n+\Pn (1) 

P — Y P — F P 

1 n+1 L n,n+l J n L n+1,71 1 n+1 

Tn+i.n and r„ jn+ i are the tunnelling rates of the island between the two states n and 
n + 1 14 . An increase in V s d leads to an increase in the tunnelling rate of the electrons at 
the drain lead. We therefore expect the tunnelling events to happen "earlier" with respect 
to the case of zero bias and the current to be more phase shifted with respect to the rf 
signal (Fig. [3]^a)). Accordingly, a maximum in phase (maximum C e ff = 0.03 fF) occurs at 
V s d = 30 /xV (FigGltc)). When V s d is increased above this value the contribution of the source 
to drain potential becomes dominant, the tunnelling events are no longer sensitive to the rf 
drive, therefore leading to a decrease in the phase shift. In Fig. E](b) and (c) we report the 
calculated phase shift as a function of both n g and V s d- The simulations have been performed 
at a finite temperature of 120mK and with dn g = 0.09e. By comparing the simulations with 
the experimental results, we observe that there is a good agreement in the n g dependence. 
The V s d dependence qualitatively agrees with the experimental observations, but presents 
some discrepancies that could be fitted by using a higher value for the temperature or dn g . 
We could not explain, however, the origin of this disagreement. 

We next consider the behavior of the amplitude. From Fig. H^b) and (c), we observe 
that the reflected (dissipated) signal reaches a minimum (maximum) at n g = 0.5 e and 
V s d = V. When V s d = V, the SET effectively behaves as a single electron box (SEB) and 
its dissipative response to a radio- frequency signal can be understood in terms of the Sisyphus 
resistance.- On average electrons tunnel after the degeneracy point between the lead and 
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FIG. 4. (Color online) (a) Reflected signal amplitude as a function of the V s d and 



rig. The decrease at V s d 



V and n g = 0.5 e is explained with an increase in 



energy dissipation by the SET due to the Sisyphus tunnelling processes. The top 
right diagrams show the different tunnelling processes involved at zero and positive 
V s d- In the first case, the Sisyphus tunnelling events lead to the absorption of extra 
energy from the resonator and thus to a decrease in the reflected power amplitude. 
In the second case, the energy needed for the tunnelling is provided by the battery 
and no decrease in amplitude is observed, (b) & (c) Graph of the measured and 
calculated variation of the reflected power amplitude as a function of the dc gate 
voltage (b) and of the source to drain voltage (c). 



island chemical potentials is passed, therefore energy dissipation occurs in each half of the 
rf-cycle, dissipating energy from the resonator. From the average power P dissipated by the 
SET (Fig. HJ), we deduce the value of the effective resistance R e ff, R e ff — ((1 / 'C ' g )dn g ) 2 /2P '. 
Our resonator is lossy (Q ~ 60), so we find P by measuring the relative variation of the 
reflected power. At V s d = V, the effective resistance reaches a minimum value R e ff = 
16MQ. 

As V s d becomes greater than the rf amplitude, the tunnelling events onto and off the island 
are driven largely by the dc bias rather than the rf excitation (Fig.|U(a))and the contribution 
from the Sisyphus dissipation becomes smaller. The absence of the CB diamonds in Fig. E](c) 
shows that the amplitude response is completely dominated by the Sisyphus dissipation, 
which overshadows the effects of the SET resistance variation at higher values of V s d- The 
average power dissipated in one period of the rf drive can be calculated as 



P = ^ / " (E n+1 - £ n )P n #[(£ n+1 - £„)P n ]cft 
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where P n +i — E n is the energy difference between n+1 and n states. The Heaviside step 
function (6) is ensures that the transitions between the states n and n + 1 are energetically 
allowed. There is close agreement between the theory and the experimentally measured 
amplitude response as a function of n g and V s a (Fig. @Jb) & (c)). 
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FIG. 5. (Color online) Calculated variation of the phase shift and of the power 
dissipated in the device, P, as a function of the ac excitation frequency fo for 



0.5e and V s d = OV. The arrow indicates the frequency at which we worked in 



our experiment. On the graph of the phase shift three points are shown: A indicates 
the region in which the electrons' tunnelling rate is comparable to the frequency 
of the rf-drive, B marks the minimum value of the phase shift, which is reached at 
higher excitation frequencies, and C indicates the region in which the frequency of 
the excitation is much higher than the electrons' tunnelling rate. The three diagrams 
on the right show the tunnelling events involved in each of these three cases on the 
time scale of the rf-drive. 

We expect that the reactive and dissipative components of the SET response to an rf- 
drive will strongly depend on the frequency of the excitation. In our experiment we don't 
have the flexibility to make a frequency dependent study of this response, however, we show 
in Fig. [5] the simulated variation of the reflected signal phase shift and of the dissipated 
power, P, as a function of the excitation frequency. As the frequency of the drive becomes 
much higher than the electrons' tunnelling rate, the tunnelling events are random and no 
longer correlated to it. However, the probability for them to occur is higher at the maximum 
amplitude of the rf-drive as the tunnelling rate is higher. Accordingly, the power dissipated 
in the device saturates towards a maximum value at higher frequencies. C e //, on the other 
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hand, tends to zero and the phase shift of the reflected signal is only determined by - 1 c ^ 9 . 
It is interesting to note that for frequencies higher than our working frequency (347 MHz), a 
negative phase shift is expected, opposite to what we observe. In diagram B. of Fig. |5j this 
is explained with the tunnelling events happening in the second half of the rf-cycle, resulting 
in a positive shift of the net current with respect to the drive. 

By connecting the gate, rather than the source or drain, to the radio-frequency circuit we 
are able to determine the Sisyphus resistance and tunnelling capacitance contributions to 
the SET impedance. These effects will also be present in a conventionally measured rf-SET 
(resonator on source-drain) and should be included if the impedance needs to be accurately 
known. 
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